Dixon's extended bodies and impulsive 
gravitational waves 



D. Bini a ' b , P. Fortini c , A. Geralico b ' d , A. Ortolan 6 

& Istituto per le Applicazioni del Calcolo "M. Picone", CNR 1-00161 Rome, Italy 

^International Center for Relativistic Astrophysics - ICR. A. 
University of Rome "la Sapienza", 1-00185 Rome, Italy 

c Department of Physics, University of Ferrara and INFN Sezione di Ferrara, 

I- 44 100 Ferrara, Italy 

d Physics Department, University of Rome "La Sapienza", 1-00185 Rome, Italy 

°INFN - National Laboratories of Legnaro, 1-35020 Legnaro (PD), Italy 

Abstract 

The "reaction" of an extended body to the passage of an exact plane gravitational 
wave is discussed following Dixon's model. The analysis performed shows several 
general features, e.g. even if initially absent, the body acquires a spin induced by 
the quadrupole structure, the center of mass moves from its initial position, as well 
as certain "spin-flip" or "spin-glitch" effects which are being observed. 
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1 Introduction 



The equations of motion for an extended body in a given gravitational back- 
ground were deduced by Dixon [1,2,3,4,5] (hereafter "Dixon's model") in mul- 
tipole approximation to any order. In the quadrupole approximation they read 



DP^ 1 



, - 'S a(S - I J^R^s^ = F* 8 ^" + F( quad ^ , (1) 

aTu 2 6 

n a 

= 2P [fl U u] + -r^R v \ a p = 2P [ »U v] + D (quad) ^ , (2) 

aTu 3 

where P M = [iU^ (with U p ■ U p — — 1) is the total four-momentum of the par- 
ticle, and is a (antisymmetric) spin tensor; U is the timelike unit tangent 
vector of the "center of mass line" C v used to make the multipole reduc- 
tion, parametrized by the proper time tjj. The tensor J 01 ^ 6 is the quadrupole 
moment of the stress-energy tensor of the body, and has the same algebraic 
symmetries as the Riemann tensor. 

In this paper we limit our considerations to Dixon's model under the fur- 
ther simplifying assumption [6,7] that the only contribution to the complete 
quadrupole moment J a ^ s stems from the mass quadrupole moment Q al3 , i.e. 
we write 

ja^s = _ w WQP]b v s\ ^ Q^U pP = . (3) 

In order the model to be mathematically consistent the following additional 
condition should be imposed [1] to the spin tensor 

S^U pu = , (4) 

to ensure the correct definition of the various multipolar terms. It is also 
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convenient to introduce the spin vector by spatial (with respect to U p ) duality 

SP = ± Va ^ s U«S, 5 , (5) 

where rja^s = V - 9 e ai3-y6 is the unit volume 4-form and e a p^ (e i23 = 1) is the 
Levi-Civita alternating symbol, as well as the scalar invariant 

s 2 = \s, u S^ = S^S" , (6) 

which is in general not constant along the path. 

Within this scheme, a general relation between U and U p can be obtained by 
using Eqs. (1), (2) and (4) (see e.g. Eq. (2.17) of Ref. [7]) 
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where 

(U-P) 

pa _|_ J_ p x JJ (quad) Act 
/i 2 /i 2 

In the case of vanishing quadrupole tensor (i.e. F^ quad ^ = = _D( c t uad )' tiy ) 
Eq. (7) reduces to the following one formerly discussed by Tod, de Felice and 
Calvani (see Eq. (14) of Ref. [8]) 

1 



ir = jU'P) 
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where 

5 = 1 4- 

4/i 



5=1 + —R^ up S^S^ . (10) 



In Dixon's model there are no evolution equations for the quadrupole moments 
of the body; therefore the system of equations (1) and (2), even if completed 
with conditions (3) and (4), is not self-consistent and one must assume that 
all unspecified quantities are known as intrinsic properties of the body under 
consideration. Moreover, the test body assumption means that mass, spin 
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and quadrupole moments (independent among each other) must all be small 
enough not to contribute significantly to the background metric. Otherwise, 
backreaction must be taken into account. 

We investigate here how a small extended body, with its center of mass initially 
at rest, reacts to the passage of an exact plane gravitational wave, extending 
previous results [9,10] limited to the spinning structure of the body. 

2 Motion of extended bodies in the spacetime of an exact gravita- 
tional plane wave 

Consider the metric of an exact plane gravitational wave propagating along 
the z direction of a coordinate frame written in the form [11] 

ds 2 = Va pdx a dxP - H(dt - dz) 2 , x a = (t, z, x, y) , (11) 

where r] a/ 3 = diag[— 1, 1, 1, 1] and 

H = hit - z)xy + hi 2 {t - z){x 2 - y 2 ) . (12) 

with h\ and h 2 two arbitrary functions associated with the two polarizations 
of the wave. Units are chosen here so that c — 1 — G. 

Let us introduce the family of static observers (with respect to the chosen coor- 
dinate system) with four- velocity = 1/ y/1 + H d t and orthonormal adapted 
spatial frame 

ej = d t + Vl + Hd z , e 2 = 9 X , e 3 = d y ■ ( 13 ) 

1 + ti 

These observers are geodesic only when located at x = = y. 
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Let the "center of mass line" be generic, i.e. with timelike unit tangent vector 
U given by 

U = 7 (e 6 + v & e & ) , 1 = (l-v 2 y l l\ 1^ = 8^, a, b =1,2, 3. 

(14) 

Analogously, let the total four-momentum P = fiU p also have the general form 
U P = 7 P (e Q + <e a ) , 7 = (1 - 4)~ 1/2 , ^ = • (15) 

Let us assume -without loss of generality- that the center of mass of the body 
is initially at rest at the origin of the coordinates, i.e. the associated world line 
Uq has parametric equations 

t = r Uo , z(r Uo ) = , x(t Uo ) = , y(r Uo ) = , (16) 

where tu denotes the proper time parameter, and hence the unit tangent 
vector reduces to Uq = e \ x=y=z= o = d t . The passage of the wave (11) will 
modify the kinematical state of the body (according to Dixon's model), due 
to the coupling between the wave and body's structure parameters. In order 
to avoid backreaction effects the natural length scales associated with the 
body, i.e. the "bare" mass /io, the natural spin length \S a \/ uq and the natural 
quadrupolar length (\Q al3 \/ Ho) 1 ^ 2 , must be small enough if compared with 
certain background scale, say Cgw (this can be associated, in turn, with the 
polarization functions h\ and hi having both the dimensions of 1/length 2 ). 
Therefore, in solving the whole set of equations (l)-(3), we will neglect terms 
which are higher order than the first in the spin as well as quadrupole length 
scales, according to 

11 = 110 + 11, u = u + u , p = hqUq + P , s^ = s^ , cr = or, 

(17) 
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the tilde denoting first order quantities which therefore must be evaluated 
along the unperturbed center of mass line (16). Hereafter, we will simply use 
the expression "to first order" to mean that spin and quadrupole terms in the 
various quantities are retained up to the first order only. 

The tangent vectors (14) and (15) thus become 

U = dt + v h d h , U p = d t + v% , (18) 

where v and v v are first order terms according to Eq. (17) and hence (7, j p ) ~ 1 
at that order. 

The conditions (4) on the coordinate components of the spin tensor imply 

S01 = S02 = S 03 = 0; (19) 

the remaining components can then be re-expressed in terms of the frame 
components (with respect to the static observer, see Eq. (13)) of the spin 
vector as 

S 23 = S\ 5 12 = - ? _, Sxs = 7 _. (20) 

Similarly, the conditions (3) 2 on the quadrupole tensor imply 

Qoo — Q01 — Q02 — Q03 — . (21) 

The spin force turns out to be 



_p(spin) _ ^(spin)2 e _ _|_ ^p(spin)3 e , 

' \hS^ - h 2 S 3 )e^ - {hS 3 + h 2 S 2 )e3 , (22) 



2(1 + //) 

to first order in the spin quantities. 
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The quadrupole force turns out to be 

^(quad) = ^uad^^uad)^. = _ { h 2 f + 2h 1 Q 23 ) fa-e^ , (23) 

to first order in the quadrupole quantities, where we have introduced the 
quantity / = Q 22 — Q 33 and the overdot means derivative with respect to 
time. 

Eqs. (f) thus reduce to the following set 

_ = ^(quad)6 ^ ji — = F (quad)i , 
dr v dr v 

drj/ dru 

so that the spin force affects only motion on the wave front, whereas the 
quadrupole force acts along the direction of propagation of the wave itself. 
The integration constants arising from the mass equation is fixed by imposing 
/2(0) = 0. The three integration constants coming from the equations for the 
linear velocities must be left indeterminate at this stage; in fact, they 
will enter the following Eq. (25) for u a , and will be fixed by requiring that 
i>«(0) = 0, according to Eq. (16). 

Consider then the evolution equations (2) for the spin tensor. By using the 
supplementary conditions (19) and Eq. (20) they give three algebraic relations 
between the spatial linear velocities u a of U and of U p 



2 ~2 i /Ol3 , u /Ol2\ 







^^^-^-(h^-hQ 12 ), (25) 
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plus three evolution equations for the spatial components S a of the spin tensor 
to be integrated together with the initial conditions S a (0) = Sfi: 



S l = --{h 1 f-h 2 Q^) , 

& = -^iQ i§ + h 2 Q^) = ~ , (26) 



Note that the overdot, used to denote derivative with respect to the coordi- 
nate time, coincides in this case with the differentiation with respect to the 
proper time r Uo , as indicated in Eqs. (16). Eqs. (26) thus imply that even if 
initially absent the spinning structure will be acquired by the body during the 
evolution, due to its quadrupolar structure. Similarly, for a purely spinning 
body the spin components are necessary constant. 

8 



Summarizing, the whole set of equation is 



(Mass) 



(Momentum) 



(CM) 



(Spin) 



ii. = -(h 2 f + 2h 1 Q 26 ) , 
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(/i 1( S 3 + h 2 S 2 ) 
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i>i + ^(h 1 Q^ + h 2 Q ii ) 



2/xq 



2/i C 
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■ (hJ-hQ 23 ) 



(hQ r2 - h 2 Q ls ) 



~(h 1 Q m + h 2 Q ii ) 



(27) 



The modification to the initial trajectory (16) of the body after the passage 
of the wave is then obtained by integrating Eqs. (27), taking into account 
that u a = dx a /drjj with initial conditions x a (0) = 0, once the profile of the 
polarization functions hi and h 2 is specified. We will explore below the two 
complementary cases of single polarization for the wave. 

2.1 Single polarization: case 1 (hi^O, h 2 = 0) 



Let us consider the case of single polarization: h 2 — 0. The whole set of 
equations (24) and (26) reduces to 
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■X _ 1 L ^i23 ~2 _ 1 . C 2 ~3 _ 1 ;. c-'i 



p 2/x p 2/i 2/i 

= - , ^ = % + , *> s = ^ + ^-h^ , 

2/i p 2/i p 2/j 

# = -^i/ , 5 s = , ^ = -^Q iS . (28) 

In order to explore the physical content of this situation we assume that all the 
nondiagonal frame components of the quadrupole tensor vanish, the remaining 
ones being constant. Such hypotheses, which can be easily released only leading 
to an increasing computational work, will allow to find out explicit solutions for 
all quantities straightforwardly, without changing the physical interpretation. 
Eqs. (28) then imply ft — and v a = and 

S l = -\hif, ^ = 0, & = | (29) 

The solution corresponding to the initial conditions /2(0) = 0, i> a = and 
S a (0) = Sq implies p, — for the mass and 

K = = v\ % = ^-H.Sl = v\ vl = —H.Sl = ifi , (30) 

for the center of mass world line, while 

S l = —HJ + Si , S* = SI, 5 s = 5j, (31) 

for the spin, where we have introduced the notation 

H 1 (r u )= Hhim- (32) 
Jo 

As a result we see that the center of mass line generally moves, but not in 
the direction of propagation of the wave. Moreover, in the case in which the 
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spin is initially absent (Sft = 0) the center of mass remains at rest, but the 
body acquires a varying spin in the direction of propagation of the wave due 
to its quadrupolar structure. It is also interesting to note that the component 
of the spin S 1 along the direction of propagation of the wave can change its 
sign if the duration of the wave is long enough, leading to a spin-flip which 
can be eventually observed. Similarly, a suitable polarization function such 
that Hiijjj) vanishes as tu — > oo guarantees that asymptotically the spin 
component S 1 goes back to its initial value, an interesting situation which will 
be sketched in the next section. 

The modified trajectory turns out to be only affected by the spinning structure 
and is given by 

t = r v , z = 0, x = — WiSg , y = --^—TixSl , (33) 

where 

H^tu) = r #!(0d£ = r d^ f hWdri . (34) 

JO JO JO 

The spatial orbit is then the line 

y = -|l . (35) 

(1) Impulsive GPW: 

In this case ^i(£) = Ai5(£), so that 

%) = y ' KM = y Tu ' (36) 

where the factor of 2 comes from integrating delta function over an in- 
terval which has as an extreme, as customary. 

(2) Sandwich wave with finite amplitude: 
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In this case hi(g) = £i[0(f) - 0(f - Co)}, so that 



B±Tu 0<T V < 7# 



(37) 



and 



\B\Tv 



< Tfr < 7# 



B^Tu - ^l(T°) 2 T U >T l 



u 



(38) 



2.2 Single polarization: case 2 (hi — 0, h 2 ^ 0) 

Let us consider the case of single polarization: h± — 0. The whole set of 
equations (24) and (26) reduces to 



4/i p 



1 1 

-?3 ~3 _ L u c2 



2/i 



h 2 S 6 , = -—h 2 S z , 

ZfiQ 



2/xo 2 /^o 2/i 

# = ~\h 2 Q™ , & = -\h 2 Q i% , S s = - l -h 2 Q li . (39) 



Let us assume again that all the nondiagonal frame components of the quadrupole 
tensor vanish, the remaining ones being constant. Eqs. (39) thus reduce to 



12 



1 . . - 1 . . - 1 - . - 1 

V = -Mf , v\ = — hf , v 2 = -—h 2 S 3 , = -—h 2 S 2 , 
4 p AfiQ 1 2/i 1 2/i 

v l = V\-^-h 2 f, & = i$, ^ = ^ p , 

Z/i 

S 1 = , ^ = , = , (40) 



with solution 



4/i 

^ = ^> ^ = 52, 5 s = 5*, (41) 



where the initial conditions /2(0) = 0, u a = and S a (0) = Sfi have been 
imposed and 

H 2 (ru)= [ TU h 2 (m. (42) 
Jo 

From the above relations we see that in this case the center of mass line 
generally moves, without specific relations with direction of propagation of 
the wave. Moreover, if the spin is initially absent (Sfi = 0), the center of mass 
moves along the direction of propagation of the wave due to the quadrupole, 
but the spin of remains zero. In addition, the body acquires a varying mass. 

The modified trajectory turns out to be 

t = W, z = ~^-H 2 f, x = ~^-H 2 Sl , y = -^-H 2 Sl , (43) 

4/i ^/io ^0 

where 

H 2 (tu) = / tf 2 (0df . (44) 
Jo 

The spatial orbit is then the line 

y = . (45) 
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The solution for this case can be found similarly to case 1 for both impulsive 
GPW (h 2 (^) = A 2 S(^)) and sandwich wave with finite amplitude (h 2 (^) = 



3 Discussion and concluding remarks 

We have studied how a small extended body, with the center of mass initially 
at rest, interacts with an incoming exact plane gravitational wave. The body 
is spinning and also endowed with a quadrupolar structure. We have discussed 
its motion by assuming that it can be described according to Dixon's model 
and by solving the corresponding set of evolution equations in the case in 
which the wave has a single polarization, for simplicity. 

A number of interesting results have been discussed. For instance, in general 
a) even if initially absent, the body acquires a spin induced by the quadrupole 
tensor; b) the center of mass moves from its initial position and the projection 
of the orbit on the wave front is a straight line, whose inclination depends 
on the initial spin of the body; c) special situations may occur in which cer- 
tain spin components change their magnitude leading to effects (e.g. spin-flip) 
which can be eventually observed. 

This interesting feature recalls the phenomenon of glitches observed in pul- 
sars: a sudden increase in the rotation frequency, often accompanied by an 
increase in slow-down rate (see e.g. [12,13,14] and references therein). Cur- 
rently, only multiple glitches of the Crab and Vela pulsars have been observed 
and studied extensively. Larger glitches in younger pulsars are usually followed 
by an exponential recovery or relaxation back toward the pre-glitch frequency, 
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while for older pulsars and small glitches the jump tends to be permanent. 
The physical mechanism triggering glitches is not well understood yet, even if 
these are commonly thought to be caused by internal processes. 

If one models a pulsar by a Dixon's extended body, then the present analy- 
sis shows that a sort of glitch can be generated by the passage of a strong 
gravitational wave, due to the pulsar quadrupole structure. In fact, from Eq. 
(31) we see that the profile of a polarization function can be suitably selected 
in order to fit observed glitches and in particular to describe the post-glitch 
behavior. 

We have just considered here a generic extended body, without exploring the 
possibility that it could actually represent a real pulsar. In fact, the inte- 
rior structure of a neutron star requires taking into account all the nuclear 
and hydrodynamical processes. This is beyond the scope of our paper. Fur- 
thermore, the observed slow-down of the period of a pulsar is expected to 
be associated with gravitational wave emission, whereas we have neglected 
backreaction effects on the background field. This analysis is enough to argue 
that the phenomena of pulsar glitches are compatible in principle with a pure 
relativistic model (Dixon's model). 

In view of getting in the next year (apparently) an enhanced phase for the in- 
terferometric detection of the gravitational waves (for both LIGO and VIRGO 
the sensitivity should be 5 times increased in comparison with the present one) 
the effects discussed here -in a very simplified form- constitute an interesting 
situation to be further explored. 
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